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We examine the prospects for detecting 7-rays from dark matter annihilation in the six most 
promising dwarf spheroidal (dSph) satellite galaxies of the Milky Way. We use recently-measured 
velocity dispersion profiles to provide a systematic investigation of the dark matter mass distribution 
of each galaxy, and show that the uncertainty in the 7-ray flux from mass modeling is less than a 
factor of ~ 5 for each dSph if we assume a smooth NFW profile. We show that Ursa Minor and Draco 
are the most promising dSphs for 7-ray detection with GLAST and other planned observatories. 
For each dSph, we investigate the flux enhancement resulting from halo substructure, and show 
that the enhancement factor relative to a smooth halo flux cannot be greater than about 100. This 
enhancement depends very weakly on the lower mass cut-off scale of the substructure mass function. 
While the amplitude of the expected flux from each dSph depends sensitively on the dark matter 
model, we show that the flux ratios between the six Sphs are known to within a factor of about 10. 
The flux ratios are also relatively insensitive to the current theoretical range of cold dark matter 
halo central slopes and substructure fractions. 

PACS numbers: 95.35.+d,14.80.Ly,98.35.Gi,98.62.Gq 



I. INTRODUCTION 

In the ACDM cosmological model, Cold Dark Matter 
(CDM) comprises approximately one- fourth of the total 
energy density of the Universe [l|]. However, the nature 
of dark matter remains unknown. Extensions to the stan- 
dard model, such as those based on supersymmetry 
and universal extra dimensions [i| , predict the existence 
of stable, weakly interacting massive particles (WIMPs) 
with mass ~ [10 1 — 10 4 ] GeV, which provide excellent can- 
didates for cold dark matter. In these models, WIMPs in- 
teract gravitationally as well as weakly, therefore WIMP 
annihilation can produce 7-ray photons. 

Present and next-generation 7-ray observatories such 
as STACEE @L HESS @, MAGIC 0], VERITAS i, 
CANGAROO H, GLAST 0, and HAWC [HI will 
search for the signatures of dark matter annihilation. 
The nearest location to search for this signal is the center 
of the Milky Way, although uncertain backgrounds from 
astrophysical sources would make the clean extraction of 
such a signal difficult [H, O E3] • Additionally, there is 
wide empirical uncertainty as to the shape of the central 
dark matter density profile, which may have been altered 
by the growth of a supermassive black hole [HI, [l6| or any 
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process which can exchange energy between the baryonic 
and dark matter components (e.g. [l7l. fl8l [T^ ) . 

In the case of dwarf spheroidal galaxies (dSphs), as- 
trophysical backgrounds and baryonic-dark matter inter- 
actions are expected to be largely absent. The Milky 
Way system contains at least 18 dSphs, which are ob- 
served to be low-luminosity systems with an extent ~ 
kpc. Based on their stellar mass to light ratios, dSphs 
contain of order O(10 1 — 10 2 ) more mass in dark matter 
than in visible light [20( and thus are ideal laboratories 
for studies that are sensitive to the distribution of dark 
matter. Furthermore, their relative proximity and high 
Galactic longitude and latitude makes them ideal for high 
signal-to-noise detection. 

In this paper, we consider the prospects for 7-ray de- 
tection from dark matter annihilation in six dSphs of 
the local group. The six dSphs are selected because of 
both their proximity and estimated masses, the latter of 
which is based on the most recent measurements of their 
velocity dispersion profiles. We estimate the range of al- 
lowable distributions of dark matter that satisfy the ob- 
served velocity dispersion profiles, and deduce the 7-ray 
flux expected from each dSph. We focus on quantifying 
the uncertainty in the predicted fluxes that comes from 
the dark matter density distribution in each system. As 
part of this uncertainty, we determine the flux contribu- 
tion of substructure within the dSph dark matter halos. 

Past work in the literature considered detecting 7-rays 
from dark matter annihilation in Milky Way-bound dark 
matter halos: dSphs were studied in fill |2U I23I [23| . 
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more massive galaxies in the local group were considered 
in [25l| potentially dark subhalos were studied in [26H27I. 
I2II l29l. l30l . l3lj ] , and the prospects of detecting microhalos 
were explored in [H, [33| • 

In comparison to previous studies of dSphs, our work is 
the first to combine theoretical predictions for CDM halo 
profile shapes and normalizations with specific dynami- 
cal constraints for each observed system. Though the 
observed velocity dispersion profiles are equally well fit 
by both central density cores and cusps, we restrict our- 
selves to inner profile shapes p cx r~ 7 with 7 ~ 0.7 — 1.2 
[34l Hil, because this is what is expected for the sub- 
set of dark matter candidates that actually annihilate 
into photons (CDM). We show that the primary uncer- 
tainty in the smooth dark matter flux contribution for 
CDM halos comes not from the relatively narrow range 
of central cusp slopes, but from the density and radius 
normalization parameters, p s and r s for the halo. As we 
show below, the published velocity dispersion data along 
with the predicted relations between p s and r s for CDM 
halos allow a tight constraint on the dark-halo density 
contribution to the annihilation signal. 

While the value of the expected flux signal for each 
dSph is sensitive to the (unknown) nature of the under- 
lying dark matter candidate, we demonstrate that the 
relative flux from system-to-system is significantly con- 
strained. Ursa Minor is the most promising dSph can- 
didate for detection and we present the expected 7-ray 
flux ratios between the remaining five dSphs and Ursa 
Minor. We also demonstrate that enhancement of the 
signal due to the presence of substructure in dSph halos 
themselves increases the predicted fluxes by at most a 
factor of - 100. 

This paper is organized as follows. In section [Til we 
discuss the 7-ray annihilation signal expected from CDM 
halos and the enhancement of the flux due to the presence 
of substructure within the dSph dark matter halos. In 
section HIll we discuss the dynamical modeling of the dSph 
galaxies. In section IIVI we present our results, and we 
conclude in section [Vl Throughout the paper, we assume 
a ACDM cosmological model with f2 m = 0.3, £!a = 0.7, 
h = 0.7 and <j 8 = 0.9. 



II. GAMMA-RAYS FROM ANNIHILATION IN 
COLD DARK MATTER HALOS 



The 7-ray flux from dark matter annihilation in a dark 
matter halo with characteristic density p s and radius r s 
at a distance T> may be written as 

f^L = ±.V[{ov),M x ,dN^/dE\ C( Ps ,r s ,V). (1) 

We have explicitly divided the flux into a term that de- 
pends only on the dark matter particle and its annihila- 
tion characteristics, V^av), M x , dN 7 /dE), and one that 
depends only on the density structure of the dark mat- 
ter halo, the distance to the halo, V, and the angular 



size over which the system is observed, C{p s , r s , T>). The 
structure quantity C is defined as 
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where the integral is performed along the line of sight over 
a solid angle Af2 = 27r(l — cos#). The term that contains 
the microscopic dark matter physics is given explicitly as 



M x 



dN 7ti (av) l 

~dE m[ 



dE. 



(3) 



Here, the mass of the dark matter particle is M x , the an- 
nihilation cross section to a final state is (crv)i, and 
the spectrum of photons emitted from dark matter anni- 
hilation to that final state is dN^^/dE. Our goal is to use 
observed velocity dispersion profiles to empirically con- 
strain the C term. This allows observations from 7— ray 
telescopes to more effectively constrain the particle na- 
ture of dark matter through V . 



A. Photon spectrum and cross sections 

As a fiducial case, we consider neutralino dark matter 
in order to determine an appropriate value for V . Neu- 
tralino annihilation to a photon final state occurs via: (1) 
loop diagrams to two photons (77), each of energy E~ n = 
M x ; (2) loop diagrams to a photon and a Z° boson (jZ°) 
with a photon energy of E lZ o — M x [l - (M z o /2M X ) 2 ]; 
and (3) through an intermediate state that subsequently 
decays and/or hadronizes, yielding photons (h). For this 
latter case, the resulting photon spectrum is a continuum 
and is well-approximated by [l2| 
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where (ai,a 2 ) = (0.73,7.76) for_WW and Z°Z° final 
states, (ai,a 2 ) = (1.0,10.7) for bb, (ai,a 2 ) = (1.1,15.1) 
for ti, and (ai,a 2 ) — (0.95,6.5) for uu. The cross sec- 
tions associated with these processes span many orders 
of magnitude. For the direct annihilation to a 77 or 7Z 
final states the maximum presently allowed value of the 
annihilation cross section to these final states is roughly 
~ (<jv) 1 ~ tl z° ~ 10~ 28 cm 3 s _1 . The total cross section 
associated with photon emission from the hadronization 
of the annihilation products has a corresponding upper 
bound of (av)h ~ 5 x 10 -26 cm 3 s~ 1 . In the most opti- 
mistic scenario, where the cross sections are fixed to their 
highest value and the mass of the neutralino is ~ 46 GeV, 
so that V = Vsvsy * 10- 28 cm 3 s- 1 GeV _2 . 

The value of V will be different for different dark mat- 
ter candidates. For example, in models of minimal uni- 
versal extra-dimensions, the annihilation cross section 
and the mass of the lightest Kaluza-Klein particle can 
be significantly higher than what we assumed here (e.g., 
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M x > 800 GeV [36|]). However, we emphasize that our 
results, which constrain the density structure of dSph's 
(and therefore C) can be rescaled to any dark matter can- 
didate that annihilates to photons, by simply multiplying 
predicted fluxes from this work with V/Vsvsy- 



B. Dark matter distribution 



Assuming a (smooth) NFW profile, the C term in 
Eq. ^ becomes 



£(p s ,r s ,V) = 2ir P y s / 
Jo 
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Dissipation-less N-body simulations show that the den- 
sity profiles of CDM halos can be characterized as 



P(r) 
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r/r s 



(5) 



where r s and p s set a radial scale and density normaliza- 
tion and 7 and 8 parameterize the inner and outer slopes 
of the distribution. For field halos, the most recent high- 
resolution simulations find 5 sa 3 works well for the outer 
slope, while 0.7 < 7 ^5 1.2 works well down to ~ 0.1% of 
halo virial radii [3J, |35| . It is currently unknown whether 
there is a "universal" 7 for every halo or if there is a scat- 
ter in 7 from halo to halo. The range quoted here char- 
acterizes the uncertainty in the theoretical prediction for 
the small-r slope, and certainly provides a conservative 
range for the halo-to-halo scatter in central slope as well. 

The structure quantity C that sets the annihilation flux 
depends primarily on the r s and p s parameters for this 
range of 7 (see discussion below) and is even less sensi- 
tive to 5. In what follows, we will fix 7 — 1 and 5 = 3 
and derive empirical constraints on the (more important) 
parameters p s and r s . Note that CDM simulations also 
predict a specific relationship between p s and r s for halos 
[e.g. H3, and this prediction is at least as robust as 
the overall shape of the profile. At the end of the next 
section, we compare our direct empirical constraints on 
the relationship between p s and r s to the expected rela- 
tionship predicted from CDM simulations and use this to 
further tighten our constraints on the dark matter struc- 
ture in the dSphs. 

With 7=1 and 6=3 the profile given in Equation [5] 
is the NFW profile and we adopt this form as the basis 
for our constraints. With the asymptotic slopes fixed, 
the values of r s and p s define the profile completely. Any 
other non-degenerate pair of halo parameters also suffice 
to characterize an NFW halo. For example, halo concen- 
tration, c = R v /r s and virial mass M, define the profile 
as well. This is a less physically relevant pair for our pur- 
poses because the virial mass is set by determining the 
extrapolated radius, R v , within which the overdensity is 
equal to the virial density, p v ~ 100/9 cr it [39j | . Given c 
and M, the value of p s is determined as p s = p v c 3 //(c) 
with /(c) = ln(l + c) — c/(l + c). A second pair of pa- 
rameters with perhaps more physical relevance is V max 
and r max - These correspond to the maximum circular 
velocity curve, V c (r) = \J GM/r, and the radius where 
the maximum occurs. Kiax is often adopted as the most 
direct characterization of the potential well depth of a 
dark matter halo, especially in the case of substructure. 



where f(0,T>,s) = \JT> 2 + s 2 — 2sT> cos8/r s , the angle 
#max defines the solid angle over which the line of sight 
(LOS) integral is performed Af2 = 2n(l — cos6* max ). In 
the particular case where a dark matter halo is at a dis- 
tance T> 3> r s , such as the case of subhalos within a dSph, 
we can rewrite Eq. ^ as 
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For an NFW profile, 90% of the flux comes within the 
region f < 1. If the angular extent of r s is less than 
the solid angle of interest, i.e., t&n~ 1 [r s /'D] < cos _1 [l — 
AD, /2tt] , Eq. jJJ reduces simply to 



£(p s ,r s ) = YP 2 s r s- 



(8) 



Typical values of the field of view of 7-ray telescopes 
are ~ 10~ 2 steradians for atmospheric Cerenkov tele- 
scopes, and ~ 2.5 steradians for space-based observato- 
ries (GLAST), with angular resolutions of ~ 1 and ~ 10 
arcminutes respectively. Note that a change in the cen- 
tral density profile slope 7 will manifest itself as a change 
in the normalization of the p 2 .r^ term Eq. ([8]). For exam- 
ple, if the inner slope is as high as 7 = 1.2, then C will be 
a factor of ~ 5.6 higher than what is stated in Eq. ((8]). 
If the profile is as shallow as 7 = 0.7, then C is smaller 
by a factor ~ 6.8 



C. Substructure and density profiles 

Dark matter halos form hierarchically, so it is expected 
that they all contain some degree of gravitationally- 
bound substructure 0, EH SI, S, El- The issue of 
dark halo substructure, or "subhalos" , is important for 
annihilation signals from dSphs for two reasons. First, 
the dark matter halos of the dSphs are "subhalos" , as 
they orbit within the virial radius of Milky Way's dark 
matter halo. We might expect this to have important 
implications for their density structure. Second, dSphs 
themselves are also expected to contain abundant sub- 
structure. This "sub-sub hierarchy" should, in principle, 
continue until we reach the low-mass cutoff scale in the 
subhalo mass function, mo ~ [10~ 13 — 1O -2 ]M , which 
is approximately set by the CDM particle free-streaming 
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scale \M IM |4J, |4g, \M [50J. In this case, we might ex- 
pect significant enhancement of the annihilation signal 
compared to the "smooth" halo assumption. 

Depending on the time of accretion and orbital evolu- 
tion, a subhalo will experience varying degrees of mass 
loss as a result of tidal interactions with the host dark 
matter halo potential. Simulations suggest that the ma- 
jority of the stripped material will be from the outer parts 
of halos. The outer slope 5 of subhalo density profiles will 
become steeper than those of field halos. However, the 
interior slope, 7, will not be altered significantly [BlT ]. 
Thus, our adopted NFW parameterization for the dSph 
dark matter density profiles is a reasonable one for de- 
termining the structure factor in the annihilation signal. 
The outer density profile slope does not affect the ex- 
pected annihilation signal and the inner slope is expected 
to remain unchanged by tidal mass loss. 

It is important to note that while the central slope 7 
is not expected to change as subhalos evolve, the nor- 
malization of the central profile does evolve, as subhalos 
monotonically lose mass, even from the central regions 
[BlT ]. One implication of this is that the relationship be- 
tween p s and r s for subhalos is altered relative to that of 
field halos. The most straightforward way to characterize 
this relationship in numerical simulations is to compare 
the Knax-^max relationship for subhalos to field halos, and 
it is found that subhalos tend to have smaller r maj , va lues 
at fixed V max such that (r s /r F ) ~ 0.7(Vs/Vf) 135 
(see also [HJ HH), where the subscripts S and F denote 
stripped and field quantities. As halos orbit within their 
parent potentials, they become less dense and their scale 
radii tend to shrink as a result of tidal interactions. We 
include this possibility when we compare our empirical 
constraints on the dSph density profiles to CDM expec- 
tations below. 



the host: BC(M) = J(dN/dm)C(m)dm, where we have 
introduced the subhalo mass function dN/dm. Unfortu- 
nately a brute-force determination of B from numerical 
simulations is not feasible at this time because the subha- 
los themselves will be filled with sub-subhalos, and this 
progression continues until the CDM cutoff scale too be- 
comes important. This requires a dynamic range of ~ 13 
orders of magnitude in halo resolution, which is far from 
the current state of the art dynamical range of numerical 
simulations. 

Our goal is to determine the expected range for B, as 
well as its dependence on too. We rely on the fact that 
subhalos tend to be less dense than halos in the field of 
the same mass. More specifically, consider the case of a 
subhalo that has experienced significant mass loss, such 
that now it has a maximum circular velocity Vmax = Vs 
that occurs at a radius r max = rg . In cases of significant 
stripping, the density profile will decline rapidly beyond 
''max [e.g. l5ll ] and the total subhalo mass will be well- 
approximated as mg ~ rgVg /G. Compare this object to 
a field halo of the same mass: M-p ~ lOpUp/G, where 
we have assumed c ~ 30 such that ~ 10% of the halo's 
virial mass is contained within r max = rp. Adopting 
the numerical simulation result quoted above, (rs/rp ) ~ 
0.7(Vs/Vf) 1 ' 35 , we can derive the relative sizes of the 
subhalo and field halo r max 's and Vmax's that give them 
the same total mass: Vs — 2.2Vf and rs ~ 2rp. At fixed 
mass we therefore expect £s/£f oc (rp/rs) 3 ~ 0.125 < 1. 

The above arguments, together with the fact that 
subhalos are expected to have less substructure than 
field halos of the same mass [I^, [55j], allow us to ob- 
tain a maximum estimate for B by conservatively as- 
suming that the total structure factor C for a subhalo 
is the same as that for a host halo of the same mass: 
^s(to) = £f(to) = C{m). Suppressing the too depen- 
dence in B, this allows us to write 



D. Substructure and flux enhancement 

Equation ((6]) assumes that the structure quantity 
C{M) in the 7— ray flux is set by a smoothly-distributed 
dark matter halo of mass M. Given the expectation for 
substructure, a more realistic formulation is that £(M) is 
set by a smooth halo component, C{M) (set by Equation 
|6]), plus a substructure component, that acts to enhance 
the flux above the smooth component expectation. It is 
useful to quantify this substructure component by intro- 
ducing a "boost" factor B: 



C{M) = [l + B(M,m )]C(M), 



(9) 



We have defined the boost such that B = is a case with 
no substructure and where all of the emission is from a 
smooth halo. The boost depends on the host dark matter 
halo mass M and, in principle, on mo, the fundamental 
subhalo cutoff scale. 

The value of B is determined by the integrated an- 
nihilation factors C{m) for subhalos of mass m within 



B{M) = 
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In the last step we have used the fact that the the sub- 
structure mass function, dN/dm, is fairly well quantified 
from N-body simulations to be a power law dN/d In to = 
A(M/m) a for m < qM, with a = 1 and q ~ 0.1 
q < 1 quantifies the fact that the subhalo mass func- 
tion cannot extend to the mass of the host itself. The 
normalization A is set by requiring a fraction / of the 
host mass M to be in subhalos with mass in the range 
qM < to < qM. Motivated by numerical simulations 
HH and semi-analytic studies [iTj, Sl[ we use / ps 0.1, 
and g w 10~ 5 to obtain A = f/ln(q/g) « 0.01 for a = 1. 

To estimate the mass dependence of C, we use the p s - 
r s relation for subhalos from the model of Bullock et 
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al. [38[ for field halos in a standard ACDM cosmol- 
ogy. As field halos are expected to be more concen- 
trated this will overestimate B(M) and is thus a con- 
servative assumption. This gives c « 33(Af/10 8 MQ)~ ' 06 
for halos of mass M < 10 8 M Q . Using the approximation 
/(c) = log(l + c) - c/(l + c) « 2.6(c/33) 04 , and that 
Ps ~ c 3 //(c), we have £ oc p 2 r z a oc A/c 2 - 2 oc Af 087 . 

We could solve for B(M) numerically with the bound- 
ary condition that B(m ) = 0. However, there is a sim- 
pler way that provides an analytic estimate for B(M). 
We note that if the upper limit of the last integral in 
Eq. (fT2|) is extended to AT, then we will have an estimate 
that will be larger that the actual B(M). Since our aim 
is to estimate how large the boost can be, this is a use- 
ful manipulation. We then differentiate Eq. (O- The 
resulting equation has an analytic solution such that we 
may write 



B(M) < A 



(M/mo)"" 7 ' 



- 1 



■ i + A 



(13) 



where we have assumed that 7' = d\n(£)/d\n(Af) is a 
constant. The sub-halo mass function will flatten off at 
smaller masses and hence this, again, is a conservative 
assumption. Note that we have imposed the boundary 
condition B(nio) — 0. This does not result in B(M) 
depending sensitively on too because a — 7' + A ~ 0.13 is 
small. For a 10 8 M© dark halo, B(M) < 41 if we choose 
too = KT 13 M , while B(M) < 2 if we choose m = 
10" 2 M m . 



III. MODELING OF DWARF SPHEROIDAL 
GALAXIES 

Twenty galaxies can be classified as residing in "subha- 
los" of the Milky Way, and eighteen of these are classified 
as dSphs. Of the eighteen dSphs, nine were discovered 
within the last two years by SDSS star counts [see e.g.[63j 
and have very low luminosities and surface brightnesses. 
We consider six of the brighter dSphs in our study: Ursa 
Minor, Draco, Sculptor, Fornax, Carina and Sextans. All 
of these galaxies have measured velocity dispersion pro- 
files based on the line-of-sight velocities of ~ 200 stars, 
which may be used to constrain their the dark matter 
halo potentials. The three remaining bright dSphs are 
Leo I, Leo II, and Sagittarius. Both Leo I and Leo II are 
too far from the Milky Way (~ 250 kpc and 205 kpc, re- 
spectively) to be detectable with 7-rays, and in the case 
of Leo II no velocity dispersion profile is published to our 
knowledge. Additionally, we do not consider Sagittarius, 
as this galaxy is known to be undergoing tidal stripping 

We assume the dSph systems to be in equilibrium and 
spherically-symmetric. Under these assumptions, the ra- 
dial component of the stellar velocity dispersion, ay, is 
linked to the total gravitational potential of the system 



via the Jeans equation, 

d{p+<7 2 ) 
dr 



-p i ,(r)V 2 (r)-2/3(r)p it a 2 



(14) 



Here p+ is the stellar density profile, the circular velocity 
is V c (r) — GM/r, and the parameter (3(r) = 1 — cr 2 /cr 2 
characterizes the difference between the radial and tan- 
gential velocity dispersions. Taking to be independent 
of radius and integrating a 2 along the line-of-sight gives 
the velocity dispersion as a function of projected radius, 

R, M 



vlos(R) 



i(R) Jr 



1-/3? 



Vr 2 - R 2 



dr. (15) 



Here, I(R) is the projected surface density of the stel- 
lar distribution, and is the three-dimensional stellar 
distribution. In Eq. (fTS"]) . cr r depends on the mass dis- 
tribution of the dark matter, and thus the parameters in 
the NFW profile p s and r s . 

The surface density of stars in all dSphs are reasonably 
well-fit by a two-component, spherically-symmetric King 
profile [66j, 



I(R) = k 



1 + — 

r 2 



-1/2 



1 r l 

r 2 



-1/2 



(16) 



where r t and r c are fitting parameters denoted as the 
tidal and core radii (see Table [TJ) , and A: is a normaliza- 
tion constant. The spherically symmetric stellar density 
can be obtained with an integral transformation of the 
surface density, 



p*(r) 



k 



7rr c [l + (r t /r c ) 2 ]3/ 2 



- cos 1 z — yl 

z 



(17) 



where z 2 = (1 + r 2 /r 2 )/(l + r 2 /r 2 ). 

Recent reductions of the photometric sensitivity in the 
extreme outer portions of dSphs show the surface den- 
sity to be falling off less sharply than expected from the 
above King profile; outside of a 'break' radius, n,, the 
surface density falls off like a power-law I(R) oc R~ 2 
57, 58, 59]. Including these variations from the King 
profile have negligible effects on the results, therefore 
for simplicity we assume the spherically-symmetric King 
profile for all dSphs. We note that for the particular case 
of Draco, recent studies have used a Plummer instead of 
a King profile, as described in fU. Using a Plummer 
profile has no effect to this calculation, because the pri- 
mary difference in the fits is in the outer regions of Draco 
where the surface density is exponentially declining. 

In order to estimate the total mass in stars and its 
contribution to the total gravitational potential, we need 
to determine the typical range of stellar mass-to-light 
ratios for dSphs. Draco was considered in Lokas et al. 
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Ursa Minor 


66 


0.29 


15 ±4 


0.30 


1.50 


15-40 


Draco 


80 


0.26 


5.5 ± 1.2 


0.18 


0.93 


15-35 


Sculptor 


79 


2.2 


8.5 ± 1.0 


0.28 


1.63 


11-19 


Fornax 


138 


15.5 


11.1 ± 2.5 


0.39 


2.71 


19-36 


Carina 


101 


0.43 


6.8 ± 1.0 


0.25 


0.86 


10-15 


Sextans 


86 


0.50 


5.8 ± 1.3 


0.40 


4.0 


6-10 



TABLE I: Properties of the dSphs used in this study. The adopted distance to each galaxy is shown in the second column. 
For reference, third and forth columns list the luminosity and central velocity disp ersion for each dwarf. The fifth and sixth 
columns give the King core and tidal radii as determined from references pA IBH! |5SL l60l . l6ll . [62| . The last column shows a 
derived result: the range of halo V ma x values that simultaneously matches the observed velocity dispersion profiles and the 
CDM theoretical normalization priors (see Fig. 2). 




0.4 0.6 0.8 
radius [kpc] 



1.2 



FIG. 1: The velocity dispersion profiles for Ursa Minor, with 
data from j57|]. The short-dashed curve shows a model with 
p s = 10 8 Mq kpc -3 and r s = 0.63 kpc, while the long-dashed 
curve depicts a model with p s = 10 7 M© kpc -3 and r s — 
3.1 kpc. Both curves have (3 = 0.6. 



[6a ] , where they quote an upper limit to the stellar mass- 
to-light ratio of ~ 3, implying a total stellar mass of 
~ 6.6 x 10 5 M Q . Though the stellar populations vary 
somewhat in all dSphs, the stellar mass-to-light ratios 
are similar [2(|. This is at the very least an order of 
magnitude below the deduced total mass in dark matter 
in all cases (see below). 

There are three empirically unconstrained parame- 
ters which determine the observed line-of-sight profile 
in Eq. ([l~5]) : 0, p s , and r s . To determine the con- 
straints on these parameters, for each dSph we construct 
a gaussian likelihood function, L oc exp(— \ 2 ), where 

X 2 = J2t( a f ~ a th,i)/ 2e i- Here °t/M is the theoreti- 
cal velocity dispersion, a l the measured dispersion in 
the I th bin, and e t is the error on <j t . The assumption 



of a gaussian likelihood function on the velocity disper- 
sion is an excellent description of the data for ~ 200 
line-of-sight velocities (69[. To construct the allowed re- 
gion, we determine the \ 2 f° r each dSph as a function 
of the three parameters /3, p s , and r s . Including 7 as 
a free parameter has minimal effect on the shape of the 
allowed region, as long as 7 is restricted in the range 
0.7— 1.2 [69(. Given L, we then integrate over the appro- 
priate range of (3 to obtain the two-dimensional likelihood 
function, L', which we use to define the likelihood ratio 
Ax 2 = —2\n{L'/L' max ). We determine the allowed re- 
gion in the p s — r s plane using A% 2 = 6.2, equivalent 
to the approximate 95% confidence level region for two 
degrees of freedom. 

Figure [1] shows an example fit for Ursa Minor, where 
we have used j3 = 0.6. The short-dashed curve has a 
maximum circular velocity, V max ~ 70kms~ 1 , and the 
long-dashed curve has V max ~ 20kms -1 . These corre- 
spond to r max <~ 0.6 kpc and ~ 20 kpc, respectively. 
This particular example highlights the degeneracy that 
currently exists with the line-of-sight velocity dispersion 
data: large V max solutions are still viable as long as they 
are accompanied by an increase in the r max . 

Figure [2] shows the allowed regions in the p s — r s plane 
for each dSph. In all of the galaxies, the minima in x 2 
is not very well-defined; there is a degeneracy along the 
axis of the allowed region. This is particularly true for 
the cases where the best-fitting value of r s occurs out- 
side the region probed by the stellar distribution. In this 
region, changes to the combination of p s — r s have very 
little impact on the dark matter distribution in the re- 
gion probed by the stars, so the allowed region actually 
extends well beyond what is shown in the Fig. [2l We 
note that if the contours are created for fixed values of /3, 
then as the value of f3 is changed, the p s and r s allowed 
region shifts along the line of degeneracy [e.g. 170T|. Thus 
our predicted C contribution changes very little whether 
we keep (3 fixed or marginalize over it (as we have done), 
especially when we demand consistency with the CDM 
model expectation for the p s -r s relation (see below). 

Though Fig. [5] shows that the combination p s — r s is 



7 




i°g 1o [p s / M k P c 



FIG. 2: The allowed region in the p 3 — r s plane for the six dSphs after marginalizing over the stellar velocity dispersion 
anisotropy parameter f3. Solid lines correspond to contours with Knax of 5, 10, 20, 40, 80, 150 kms . Long-dashed lines 
represent the p 3 — r s relation as derived from the field halo relation, and the 2-er scatter above the median concentration vs. 
mass relation. Dot-dashed lines represent the p s — r s relation as derived from the tidally-stripped halo relation, and the 2-a 
scatter below the median concentration vs. mass relation. 



not well-defined, in all of the cases the data does approx- 
imately fix the density at the mean radii r* of the stellar 
distribution [6!|. Calculating the total mass of the dark 
matter within this characteristic radius, for all galaxies 
the minimum implied dark matter mass ~ 10 7 M Q , which 
occurs for the lowest implied values of p s — r s in each case. 
This is at least an order of magnitude greater than the 
contribution to the total mass in stars in all cases. 

Over-plotted in Fig. [2] are lines of constant V max in 
the p s — r s plane. Phrasing the dark matter halo prop- 
erties in terms of Kn ax allows for a direct comparison 
to CDM models, which provide predictions for the cu- 
mulative number distribution of halos at a given V mSLX . 
Although the high V max solutions are plausible by consid- 
ering the data alone, comparison to CDM models show 
that it is improbable that all of these halos have Vmax 
in the high end of the allowed regime [S H [zl Q 
(although this solution may be viable for some smaller 
fraction [lH). Typical CDM halos have ~ 1 system as 
large as ~ 60 km s . 

Dashed and dash-dotted lines in Fig. [5] enclose the 
predicted p s -r s relation (including scatter) for cold dark 
matter halos as determined from numerical simulations. 
In order to provide a conservative range for the CDM ex- 
pectation, the upper (long-dashed) lines are obtained us- 



ing the relation that is 2-a above the median for field ha- 
los in ACDM [38| and the lower (dash-dotted) lines show 
the relation implied by the the tidally-stripped Knax-^max 
relation with a 2-a scatter below the median c(M) rela- 
tion [53|]. We consider both the field and stripped 
relation because the degree of tidal stripping experienced 
by each dSph is uncertain, depending sensitively on the 
precise orbital information and/or redshift of accretion, 
two quantities that set the amount of tidal mass loss . 
The region where the CDM predictions cross with the 
observationally-allowed values of p s and r s in Fig. [5] de- 
fines a preferred model for the structure of these dark 
matter halos within the context of CDM. 



IV. FLUXES FROM DWARF SPHEROIDAL 
GALAXIES 

A. Smooth Halo 

The flux of 7-rays originating from the annihilation 
of dark matter particles is sensitive to p 2 s r 3 s (recall that 
C r-j p 2 s r^, see also Eq. (JBJ)). Even though p s and r s indi- 
vidually can vary by orders of magnitude and still satisfy 
the observed velocity dispersion profile (see Fig. [5]), the 
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FIG. 3: The allowed region in the p s r 3 — p s plane for each dSph after marginalizing over the stellar velocity dispersion anisotropy 
parameter (3 filled region, and contour levels for the expected 7-ray flux. Contours are shown for log 10 [diV 7 /cL4dt] = —13, —12, 
— 11, & —10, where the flux is measured in photons cm _2 s _1 . Solid contours depict the flux expected within a region of radius 
2 degrees centered on the dwarf, while dot-dashed contours depict the same flux thresholds for a region of radius 0.1 degree 
centered on the dwarf. The hatched regions represent the preferred region from CDM theoretical modeling (see Fig. [2] and 
discussion in text). 



product p 2 r^ is tightly constrained. Fig. [3] shows the al- 
lowed region in the p 2 r^ — p s plane. The tight constraint 
makes the predictions for 7-ray fluxes more robust. The 
hatched regions correspond to solutions that overlap with 
the CDM expectation in this parameter space. With the 
CDM prior imposed, the p 2 r^ s quantity is constrained to 
within a factor of ~ 3 — 6 in all cases. This corresponds 
to the width of the hatched regions in Fig. [3l 

The mild change in p 2 r^ with p s can be explained by 
looking at the radial velocity dispersion measure. First 
note that the contours in Fig. [3] have a common shape. 
They slope gently down for about a decade in p s , re- 
main constant and then start to slope upwards (as p s is 
increased). In addition the area where the contours re- 
main roughly constant is where r s ~ r t ; recall that r t is 
the tidal radius of the stars. The LOS velocity disper- 
sion is a weighted average of the radial velocity disper- 
sion for r < r t , see Eq. 1151 Thus the p s -r s scaling must 
trace back to the scaling of cr 2 {r). For a constant stellar 
anisotropy, we have 

fJ r -20 roc 

4(r) = — 7-T- / dr>p4r>y 2 ^ 2 M(r>) , (18) 
where M(r) is the total mass profile, which is to a good 



approximation the dark matter mass profile. We are in- 
terested in the scaling of a 2 (r) with p s and r s . It is clear 
from Eq. (fT5)l that a 2 (r) cx p s always. To understand 
the scaling with r s , we consider the following integral, 
r 3 s f r l dr' p+(r')r' 2 P- 2 (ln(l + r'/r s ) - r'/(r' + r s )). This 
integral, obtained from Eq. (fTBj) . has all the information 
about the scaling of <J 2 (r) with r s . For r s ^> r t , the 
NFW mass term can be replaced with r' 2 r Sl which im- 
plies that <J 2 {r) scales as p s r s for a given r, and hence 
we would predict that p s scales as 1 jr s for r s 3> r t . This 
prediction is verified by the shape of the contours for the 
smaller values of p s (larger values of r s ) in Figs. [21 EJ 

For r s <C r t , the NFW mass term varies slowly com- 
pared to p+ and possibly the (3 terms. This means that 
of (r) should scale (in the above limit) as p s r^ or the 
contour of p 2 r^. should increase linearly with p s . We can 
discern this behavior for Sextans, which has a large rt- 
In the intermediate region where r s ~ r t , o 2 (r) varies 
faster than linearly with r s . To see this, we first note 
that 8{r) = d In M{r)/d In r is 1.3 at r — r s and is 1.8 
at r = 0.2r s . Also, r c /r t is between 0.1 and 0.3 for 
these galaxies and most of the contribution to the LOS 
dispersion comes from the region around r c . Therefore 
we expect the contours to be along curves of constant 
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FIG. 4: Examples of the flux spectrum of Ursa Minor for 
three cases where the quanitities (log 10 p s , log 10 r 3 , M x ) take 
the values of (7.4,0.033,46) depicted with the long-dashed 
line, (7.9,-0.067,46) shown as a short-dashed line, and (7.9,- 
0.067,500) shown as the dot-dashed line. The value of C 
that corresponds to these 3 cases is [2.08 x 10 14 ,1.25 x 
10 15 , 1.25 x 10 15 ]GeVcm" 2 s _1 respectively. The units for p s 
are M0kpc -3 , while r B is in kpc and M x in GeV. No en- 
hancement of flux from substructure is included; substructure 
could increase the flux by up to a factor of 100, increasing the 
prospects for detection. The calculated flux is integrated over 
an angular region of radius 0.1 degrees centered on the dSph, 
and the value of V = Vsvsy ~ 10" 28 cm 3 s _1 GeV" 2 , which 
corresponds to the most optimistic scenario for supersymmet- 
ric dark matter (see Sec. HIJ. Open squares show the am- 
plitude of the 7— ray extragalactic emission [74|, while filled 
circles correspond to the galactic emission of 7-rays at high 
galactic latitudes [75| . 



Psf's S<J ~ C \ which explains the flat parts of the contours 
in Fig. [3J (where they overlap with the CDM priors). 

The contours in p 2 r 3 — p s plane in Fig. [3J have been 
calculated assuming a smooth dark matter distribution. 
The flux enhancement due to the presence of substruc- 
ture will be discussed below. We have used V = Vsvsy — 
10 _28 cm 3 s _1 GcV -2 , which corresponds to the most opti- 
mistic scenario for neutralino CDM. For other dark mat- 
ter candidates, the fluxes shown should be rescaled by a 
factor of V/Vsusy- 

Figure [3J shows that the most promising candidates for 
detection are Ursa Minor & Draco, with the largest flux 
coming from Ursa Minor. These two dSphs have fluxes 
~ 1CP 11 cm s , within potential reach of upcoming 7- 
ray detectors. For example, the integral sensitivity for 
a 5-cj detection in 5 years of exposure with GLAST in 
the signal dominated regime (energies above ~ 5 GeV) 
is ~ 3 x 10 _11 cm _2 s~ 1 , and therefore these two dSphs 
should be prime targets for observation with GLAST. 



The various lines in Fig. [3J show flux levels for different 
solid angles of integration centered on the dSph. Because 
most of the flux from a dark matter halo described with 
an NFW profile originates from the region inside of r s , 
integrating over an area that is larger than the appar- 
ent angular extent subtended by r s does not lead to a 
marginal increase in the flux (see e.g. Eq. ([5])). For a 
dSph at a distance T> this angular extent is tan -1 [r s /T>]. 
Integrating over an angular area which has a apparent 
radius smaller than r s leads to a reduction in flux (see 
e.g. Eq. ([7])). This is shown with the dot-dashed contours 
in Fig. [3J where the solid angle is 0.1 degrees relative to 
the solid contours which are for 2 degrees. 

In order to quantify the prospects for detection we con- 
sider the following examples. If a region of radius 0.1 
degrees centered on Ursa Minor is integrated upon with 
GLAST (with an orbit-averaged area of A c s ~ 2x 10 3 cm 2 
[Z|) for 5 years, and V » 10" 28 cm 3 s- 1 GeV -2 , then the 
range in the number of photons expected is ~ [5 — 35] 
based on the allowed range of values in the p s — r s plane. 
Integrating over the same region with a Cerenkov de- 
tector (such as VERITAS (atmospheric) or HAWC (wa- 
ter)) has the advantage of a much larger effective area 
(A e s ~ 10 8 cm 2 ), but the disadvantage of a much larger 
background (due the hadronization of cosmic rays) and 
much smaller integration timescale (of order hours in- 
stead of years). For ground detectors such as VERITAS, 
or HAWC, with an effective area .4 e ff ~ 10 8 cm 2 , and as 
an example, 50 hours of integration, the corresponding 
range in the number of photons expected is [10-70]. For 
this latter estimate we assume V ~ 10 _31 cm 3 s _1 GeV -2 
which corresponds to a neutralino of M x ~ 200 GeV and 
a threshold energy of ~ 100 GeV. 

As can be inferred from Fig. [3J the predicted fluxes 
are roughly similar to within two orders of magnitude. 
If a 7-ray flux is detected in the direction of, for exam- 
ple, Ursa Minor, then from the allowed p s — r s parameter 
space of Ursa Minor we can determine the range of ex- 
pected fluxes for the remaining five dSphs, by taking into 
account the respective allowed p s — r s parameter space 
in each case. Table Ql] provides the flux ratios expected 
relative to a flux measurement from Ursa Minor. We cal- 
culate these flux ratios by considering the highest and 
lowest flux in the p 2 r 3 — p s parameter space which is also 
consistent with the CDM priors (shaded areas in Fig.[3J. 
If the highest flux predicted from Ursa Minor is $umi 
and the minimum is $umi> t nen the- range of flux ratios 
from the rest dSphs relative to the flux from Ursa Minor 
is ^ph/^UMi " *dSph/ $ UMi- We calculate flux ratios 
for two different angular integrations, such that combi- 
nations of the two removes any correlations between the 
allowed regions by the inclusions of the distance to each 
dwarf, i.e., a same allowed value of p 2 r 3 in for example 
two different dSphs does not necessarily correspond to 
the same flux (recall that the angular extent of r s for a 
dSph at D is tan - 1 [r s /D], and that the flux is propor- 
tional to plr^g). This prediction is quite robust. First, 
because measurement of 7-ray fluxes must fall within 
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FIG. 5: Left: The predicted substructure boost factors assuming a subhalo mass function scaling as dN/dM ~ M 19 . Right: 
The dependence of the overall substructure boost factor to the slope of the subhalo mass function for Ursa Minor. The assumed 
parameters are M = 3.02 x 1O 8 M , p a = 6.3 x 10 7 M Q kpc" 3 , and r 3 = 0.8 kpc. 



the prescribed tight range in the p 2 s r^ s — p s plane if all 
dSphs are composed of the same CDM particles, and 
second because the predicted flux ratios are less sensi- 
tive to astrophysical processes, which may contaminate 
the 7-ray emission. Therefore, correlated fluxes between 
dSphs should be expected in future measurements. 

The results presented in Fig. [3J are fluxes integrated 
over an energy regime that contains the signal from dark 
matter annihilation for a fiducial value of V . Extraction 
of a signal from dark matter annihilation will depend on 
the shape of both the input and the background spectra. 

In Fig. [4J we show an example of the spectrum from 
dark matter annihilations in Ursa Minor, 
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= 5 x 10" 



dE 



5 x lO^cri^s -1 
46 GeV 
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~B + l 




1 



C(p s ,r s ,V) 



1.25 x 10 15 GeV^cm" 5 



(19) 



The assumed values of p s and r s in this example are 
consistent with the observed velocity dispersion-derived 
profile, and also with the predictions of CDM theory (see 
Fig. [2]) • This highlights the range of values that the am- 
plitude of the spectrum may take. As shown in Fig. [3J 
the structure quantity C can be a factor of 3 greater, or a 
factor of 2 smaller than what is assumed in this example. 
On the other hand, the cross section assumed represents 
an upper bound for the continuum emission of 7-rays 
from the hadronization of the annihilation products, and 
can be smaller by up to 6 orders of magnitude, while the 



dSph 



within 0.1 deg 



within 2 deg 



Draco 
Sculptor 
Fornax 
Carina 
Sextans 



0.1-3.2 
0.07-1.6 
0.07-2.2 
0.04-1.0 
0.02-0.5 



0.1-2.8 
0.05-0.7 
0.05-1.1 
0.02-0.4 
0.007-0.02 



TABLE II: The predicted flux ratios for dSphs relative to the 
7-ray flux from Ursa Minor in CDM theory. 



mass of the dark matter particle represents the lower ex- 
perimental bound for neutralino dark matter. Note that 
the boost factor B in this particular example is taken to 
be 0. We will discuss this in the next section. 

In Fig. [4] we show the spectrum of photon emission 
from within an angular region of 0.1 degrees for different 
choices of the values of p s , r s and M x . Shown are also 
the 7-ray background spectrum that have extragalactic 
[74| and galactic [75[ origin. Increasing the angular ac- 
ceptance of a detector from 0.1 degrees to radii larger 
than the projected angular size of r s does not lead to 
a significant increase in the flux from the dSph, but it 
does increase considerably the flux from the two diffuse 
components. This can be understood in the following 
way. For an NFW profile, the majority of flux originates 
from the region within r s . The distance to Ursa Minor is 
T> = 66 kpc, and therefore the observed angular extent 
of the scale radius of Ursa Minor is t&TT 1 [r s /T)] ps 0.7 
degrees. Thus, integrating over a region greater than 0.7 
degrees does not lead a substantial increase in the flux. 
On the other hand, we can find the decrease in the mea- 
sured flux that results from an integration over an area 
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dSph 


Radius of the area of 
90% flux emission in degrees 


Ursa Minor 


0.4-2.7 


Draco 


0.3-1.8 


Sculptor 


0.2-0.9 


Fornax 


0.2-1.0 


Carina 


0.1-0.6 


Sextans 


0.1-0.4 



TABLE III: The CDM-predicted angular extent in degrees 
where at least 90% of the 7-ray flux should originate for each 
dSph. 



smaller than the angular extent of r s , say, 0.1 degrees. In 
this case, we use Eq. ([7]), and find that the emitted flux 
should be a factor of 1.7 less than the value obtained 
by integrating out to 0.7 degrees. 

We therefore emphasize that increasing the area of in- 
tegration does not significantly increase the flux from the 
dSph as long as it is of order or larger than tan _1 [r s /2?]. 
However, increasing the integration area does increase 
the photon counts that originate from contaminating 
sources. In Table IIIII we show the range of the angular 
extent in degrees where 90% of the flux will originate for 
each dSph. We calculate this quantity using the values 
of r s consistent with CDM predictions (see Fig. [J). Fu- 
ture observations of the six dSphs presented here should 
concentrate on integrating over areas with radii as shown 
in Table Hill centered on each dSph. 



B. Including Substructure 



(see Fig. [2]). We note two important results: 1) the value 
of B is only weakly dependent on the cut-off scale of the 
subhalo mass function, and 2) the boost factor can take 
values which are at most of order ~ 100 in all dSphs. We 
can understand these effects by recalling the solution to 
Eq. (fl3|) . where the boost factor of a halo of mass M (with 
a subhalo mass function with a cut-off at a scale mo) is 
approximately given by B(M, mo) ~ 0.1((M/mo) ' 13 — 
1). The weak dependence on mo is a result of the flatness 
of the relationship between the concentration and mass 
in CDM halos (which is an outcome of the flatness of the 
dark matter power spectrum). For scales ~ 10 7 M Q , the 
concentration parameter scales with mass asc~ m -0 06 , 
while for much smaller scales, e.g. ~ 10 -5 M Q , it becomes 
even shallower, c ~ m~ 037 . As such, the boost factor in 
dark matter halos does not increase dramatically when 
the cut-off in the subhalo mass function is reduced. 

In the right panel of Fig. [5] we show the effect of 
the subhalo mass function power law to the boost fac- 
tor. As an example, we use Ursa Minor, with a mass 
of M = 3.02 x 10 8 M©, and characteristic density and 
radius of p s = 6.3 x 10 7 M Q kpc -3 and r s — 0.8 kpc re- 
spectively. A change in the power law index leads to 
significantly different behavior. For example, if the sub- 
halo mass function has a cut-off at 10~ 5 M©, then a 5% 
uncertainty in the subhalo mass function power law man- 
ifests itself into a difference in B by as much as a factor 
of 80. In addition, note that for a subhalo mass function 
dN/d\nm ~ m~ a , and a luminosity of C ~ p 2 s r^ s ~ m ' 87 , 
the luminosity per logarithmic mass interval in substruc- 
ture is dL/d\nm ~ £(m) dN/dlnm ~ M°- 87 ~ a . There- 
fore, for mass functions with a ~ 0.9, the contribution to 
the boost factor per logarithmic mass interval is a very 
weak function of subhalo mass. 



In this section, we study the effect of substructure on 
the 7-ray flux. We assume the same constraints on the 
dark halos as shown in Fig. O and combine them with 
Eq. (|12p in order to determine the substructure flux boost 
factor B introduced in Eq. ([9|). We will first need an 
estimate for the total mass of the dark halos for a given 
point in the p s — r s parameter space. For these estimates, 
we first use the Jacobi approximation to determine the 
tidal radius, r t ~ V(m(r t )/3/M M w(P )) 1 / 3 . Here V is 
the distance to the dSph, Mmw(1^) is the extrapolated 
mass of the Milky Way at that distance, and m is the 
mass of the dSph which we wish to determine. For the 
Milky Way mass, we use an asymptotically flat rotation 
curve, with U ma x — 220 kms~ l . Because T> is much larger 
than the typical scale radius of the Milky Way halo when 
fit by an NFW profile, our results are insensitive to the 
choice of the mass model for the Milky Way. For each 
galaxy, we thus solve for r t , and then the total mass 
within r t , given an input pair of p s — r s . 

In the left panel of Fig. [5] we show the range of values 
of the substructure boost factor for each dSph, based 
on the allowed region from combining CDM theory and 
velocity dispersion data in the p s — r s parameter space 



V. SUMMARY 

We address the prospects for detecting dark matter 
annihilation from six dwarf spheroidal satellites of the 
Milky Way. Using the stellar velocity dispersion profiles 
for each dSph, and assuming an NFW profile for the dark 
matter, we deduce constraints on both the characteristic 
density, p s , and characteristic radius, r s . We show that 
each dSph exhibits a degeneracy in the p s — r s param- 
eter space. We have assumed that the stellar velocity 
dispersion has a constant anisotropy and allowed it to 
vary. However, the degeneracy exists even if the anisop- 
tropy is kept fixed. The p s — r s degeneracy translates to a 
degeneracy in the more observationally-relevant parame- 
ters of U max and the radius where the maximum rotation 
speed is attained r max . The degeneracy direction is such 
that larger values of U max are allowed as long as they 
are accompanied by the corresponding increase in r max . 
However, this degeneracy is broken substantially because 
in CDM theory, there is a relation between p s and r s . We 
find that imposing this CDM "prior" constrains 5 km s _1 
< Kiiax < 40 km s _1 in all the dSphs we consider (see 
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Table 2). 

Assuming a smooth dark matter distribution in the 
dSph halos, we find that Ursa Minor and Draco are the 
most promising dSph's for detecting products of dark 
matter annihilation. Fornax and Sculptor are a factor 
of ~ 10 fainter, while Carina and Sextans are fainter by 
a factor of ~ 100. In the most optimistic scenario for 
neutralino dark matter, the largest-predicted flux from 
Ursa Minor is ~ 3 x 10~ n cm -2 s _1 . This is the flux 
within a 0.1 degree radius centered on Ursa Minor, and is 
within the sensitivity threshold of future detectors, such 
as GLAST [nj. Given the fact that all dSphs will have 
the same spectrum from dark matter annihilation, the 
prospects of detection may be further enhanced by stack- 
ing the signal from all 6 dSphs presented in this work. 
This can lead to an increase in the total flux up to a fac- 
tor of 2. The flux predictions presented here can easily 
be rescaled to any dark matter candidate that annihilates 
to photons (see Sec. IIVI) . 

The dark matter distribution is certainly not smooth 
and the presence of rich substructure in dark matter ha- 
los can enhance the flux from annihilation of dark matter 
particles. We calculate this enhancement resulting from 
substructure, and find that it can be at most ~ 100, inde- 
pendent of the cut-off scale in the subhalo mass function. 
In the most optimistic particle physics scenario, this en- 
hancement puts the fluxes from all the dSph's we consider 
above the threshold of future 7-ray detectors. 

While the allowed region in p s — r s parameter space is 
degenerate, we show that the corresponding range in the 
product p 2 s r^ is much more tightly constrained. This is 
important because the 7-ray luminosity from dark matter 
annihilation is L ~ Pp^, implying that the range of pre- 
dicted fluxes is narrow. We find that the observationally 
deduced values in the fluxes can vary by a factor of ~ 10 
and imposing the CDM prior further reduces the uncer- 



tainty to a factor of ~ 3 — 6, depending on the particular 
distribution of dark matter in each dSph. This range will 
only be reduced with the inclusion of more stars in the 
analysis of the line-of-sight velocity dispersion profiles. 

Throughout this work, we have assumed that the dark 
matter density profiles are described with an inner slope 
of 7 = 1 (NFW). Dark matter annihilation signal is pro- 
portional to the square of the density, so the predicted 
flux is sensitive to the value of 7. Varying the inner slope 
within the current theoretical uncertainty 0.7 < 7 < 1.2 
results in a flux increase or decrease by a factor of ~ 6. 
It is not yet clear if the spread in 7 we quote above is 
truly the scatter from halo to halo or if much of it reflects 
numerical resolution issues. If there is a distribution of 
values in 7 as large as that quoted above, and it is in- 
dependent of host halo mass, then this uncertainty will 
have to be factored into the flux predictions. If the inner 
slope correlates with mass or if the true scatter in 7 from 
halo to halo is small, then our predictions for flux ratios 
are robust. Future N-body simulations will be crucial in 
constraining the theoretical uncertainty in 7. 
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